= 
© 
© 
N 

" 

O 
© 
e 
E 
C 
© 
— 
= 


B-meson dileptonic decays enhanced by supersymmetry with large tan 8 


Zhaohua Xiong %®° and Jin Min Yang ^ 
^ CCAST (World Laboratory), P.O.Box 8730, Beijing 100080, China 
* Institute of Theoretical Physics, Academia Sinica, Beijing 100080, China 
° Institute of High Energy Physics, Academia Sinica, Beijing 100039, China 
(December 22, 2013) 


We examined the rare decays B — X,£*6(^ and B, — 4+4 y in the minimal supersymmetric 
model with large tan 8. Taking into account the gluino-loop and neutralino-loop effects, we found 
that for a large tan 8 the neutral Higgs exchanging diagrams could enhance Br(B — X,r*7 ) by a 
factor of 5 and Br(B, — rry) by a couple of orders in some part of supersymmetric parameter 
space allowed by current experiments such as b — sy, B — KMete- and B, > £*£-. The 
forward-backward asymmetry and the distributions of differential branching ratios are also found to 
differ significantly from the standard model results. Such enhanced branching ratios reach the level 
of 10^? and thus might be observable in the new generation of B experiments. 
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I. INTRODUCTION 


Flavor-changing neutral-currents (FCNC) induced B-meson rare decays provide an ideal opportunity for extracting 
information about the fundamental parameters of the standard model (SM), testing the SM predictions at loop level 
and probing possible new physics. After the observation of the penguin-induced decay B — X, and the corresponding 
exclusive channels such as B —^ K*y m rare B-decays have begun to play an important role in the phenomenology 
of particle physics. The latest measured decay ratio for B — X,y by CLEO and BELLE B is in good agreement 
with the SM prediction, putting strong constraints on its various extensions and therefore stimulating the study of 
radiative rare B-meson decays with a new momentum. 

Among rare B-meson decays, B, — Lly (€ =e, p, T) are of special interest due to their relative cleanliness and 
sensitivity to models beyond the SM BH. Since in these processes a photon is emitted in additional to the lepton 
pair, no helicity suppression exists and “large” branching ratio is expected. Other interesting decay modes in this 
context are the inclusive transitions B —^ X,£*4^. Although these rare decays have not been observed, their detection 
is expected at the B-factories which are currently running. 

These decays have been studied in the SM ra and recently, to reduce the theoretical uncertainties, the next-to-next- 
to leading order(NNLO) corrections were completed m. New physics effects in these decays have also been studied 
in some models, such as the minimal supersymmetric model (MSSM) BH, the two Higgs doublet model(2HDM) 
[L2HL5] and the technicolor model Bl. 

It is noticeable that in the SM the matrix elements of B — X,4*(^ and B, — €+£~7 are strongly suppressed by 
a factor me/my and the contributions from exchanging neutral Higgs boson can be safely neglected. In the MSSM 
im the situation is different, specially in the case of / = 7 with large tan. In this model, as studied in [m the 
contributions from exchanging neutral Higgs bosons are enhanced roughly by a factor tan? 8 and may no longer be 
negligible since currently a large tan 9 is favored both by LEP experiments and by the supersymmetry (SUSY) 
explanation fis} of the measurement for the muon anomalous magnetic moment [19]. 

We note that the most previous studies on the contributions from exchanging neutral Higgs bosons [911] mainly 
focused on the charged-current loop effects; the effects of neutral-current loops (NCL), such as gluino-loop and 
neutralino-loop, have been considered only for BY, — u"a using mass insert approximation method in [B]. A 
detail general calculation for such NCL effects is necessary. The NCL can be induced via the flavor mixing of down- 
type squarks and might be important for the following reasons. Firstly, such flavor mixings of sfermions are almost 
unavoidable in supersymmetric models. In fact, in the framework of MSSM sfermions may have arbitrary flavor 
mixings in the soft breaking terms; while in some constrained MSSM, such as low-energy supergravity models, the 
flavor mixings at weak scale could be naturally generated through renormalization equation even the flavor diagonality 
is assumed at the Planck scale o. Secondly, the flavor mixings between the third and the second generation squarks 
are subject to no strong low-energy constraints like K} — K? mixing. Thirdly, the large tan B will give rise to large 
mass splitting between two mass eigenstates of sbottoms, making the lighter sbottom (b1) even lighter. 

In this article, we will present a complete calculation of MSSM effects in the decays B, —^ (* £^» and B —^ X,U*£7 
(€ =e, u, T), taking into account the contributions from the neutral Higgs exchange with NCL. We will evaluate 


the effects on branching ratios, the forward-backward asymmetry as well as the distributions of differential branching 
ratios. In Section in, we will give a brief description of the squark mixing in the MSSM. The detailed calculations 
for Wilson coefficients of scalar and pseudo-scalar operators are included in Section mi The decays B > X,£*£^ 
and B, — £*£^^ are calculated in Section a Experimental constraints on the parameter space of the MSSM are 
discussed in Section M Some numerical results are presented in Section vil Finally, in Section [vri we give our 
conclusion. 


II. SQUARK MIXINGS IN MSSM 


The squark mass terms arise from the scalar potential which contains the supersymmetric F-term and D-term as 
well as the soft SUSY breaking terms. These soft breaking terms may have arbitrary flavor mixings. As a result, the 
squark mass terms in flavor basis, i.e., U = (ür,Cr,tr,üg, Cg, tg) and D = (dr, 8r, br, dg, Sn, bg), take the forms 


um - Mọ + mz, + m3 (3 — ey sin? Oy) cos 28 mu (Au — u cot B) (2.1) 
US mu (Ay — p cot 3) Mz + m2, +m ey sin? 0w cos 23 |’ f 
P M$ m? — mZ(3 + ep sin? Oy) cos 28 mp(Ay — utan f) (2.2) 
b mp(Aq — utan 8) M2 t mj t mZep sin? Ow cos2 |’ ; 
where ey = 2/3, ep = —1/3, my(p) is the 3 x 3 diagonal mass matrix for up(down)-type quarks. M2, M? and 
U(D) ong 


M are soft-breaking mass terms for left-handed squark doublet Q, right-handed up and down squarks, respectively. 
Au (Aq) is the coefficient of the trilinear term H2QU (H,QD) in soft-breaking terms and tan  — v/v, is ratio of 
the vacuum expectation values of the two Higgs doublets. The hermitian matrices M? p can be diagonalized by 
the unitary relations, which transfer the interaction eigenstates into the physical mass eigenstates. So in the general 
MSSM, without knowing the mechanism of SUSY breaking, squarks could have arbitrary flavor mixings. 

However, the flavor mixings in the first two generations are subject to strong phenomenological constraints, such 
as K? — K? mixing. So we only consider the flavor mixings between the second and third generations, i.e., between b 
and 8. Further, like the analysis in p3, we suppose the tree level Lagrangian is flavor diagonal and the flavor mixing 
is induced via loops. The dominant effects are from the logarithmic divergences caused by soft breaking terms. Such 
divergences must be subtracted using a soft counter-term at the SUSY breaking scale, such as Planck scale Mp. Thus 
a large logarithm factor In( M? /mi,) ~ 80 remains after renormalization. In the approximation of neglecting the 
strange quark mass, $g does not mix with sbottoms. The mixing of sz with sbottoms results in the physical states 
is given by 


81 SL 1 e eh SL 
bi )=T? | bs |=| a cosh sind br |, (2.3) 
bs bn €? —sin@, cos bn 
where 
pa (uà + M2+M2 + Ad?) cos Øs + ms A? sin 0s 
€1 = —6b—35 ee ——————— — — —ÁÀdJa4 
sinf 8 mi ms, 
me [7 (uà - MEME |4412) sin 0, + m, A; cos 05 
2 — —Cb— — Ls 
iia tite 
€; = —€ cos Oy + ea sin Mp, 
€) = —€ sin Oy — €2 cos Ob. (2.4) 
: : " ‘ K* Ky M2 
Here 6, is the left-right sbottom mixing angle, K;j are the CKM-matrix elements, and c; = = can tap Oe 
WwW 


Note the mixing matrix TU between the left-handed scharm and stops takes the similar form as T?, and under the 
assumption that the flavor mixing between the second and third generation squarks is at least one order lower than 
the third left- and right-hand squarks mixing, the rotation matrix TP are approximately unitarity. 


III. CALCULATION OF SCALAR AND PSEUDO-SCALAR WILSON COEFFICIENTS 


In the MSSM, the short distance contribution to b — s£* £^ decay can be computed in the framework of the QCD 
corrected effective weak Hamiltonian, obtained by integrating out heavy particles, i.e., top quark, WF, Z bosons in 
the SM and the sparticles, 


10 
Hops = “SE Rais Yo (Ci) Oslin) + Co, (ue) Qiu) (3.1) 


i=1 


where O;, Q; are operators given in Eqs. (7.1) IE (7.4) [L3] and C;, Co, are Wilson coefficients renormalized at the 
scale ur t] 
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FIG. 1. The Feynman diagrams involving charged Higgs which give dominate contribution to Ca, 2 


In the MSSM the additional contributions to operators in Eq. (B.1) can be characterized by the values of the 
coefficients C; and Cg, at the perturbative scale my. For the processes we will study, it is only relevant with the 
effective Wilson coefficients C7,9,19 which have been computed in Ref. and Co, , of additional scalar and pseudo- 
scalar operators. In this section we will focus our attention on the calculation of the Wilson coefficients CQ, , (mw) 
with the assumption that except for the third generation squarks all sfermions are degenerate and have a mass of ~ 1 
TeV. Besides the contribution from box diagrams, the neutral Higgs-bosons exchange diagrams include totally five 
classes of loops: (1) W boson and up-type quarks, (2) charged Higgs and up-type quarks, (3) charginos and up-type 
squarks, (4) neutralinos and down-type squarks, and (5) gluinos and down-type squarks. As pointed out in Section Il 
the last two classes of loops have not been calculated in the literatures. Now we take all of them into account and 
use the Feynman rules presented in [io]. Since we are only interested in large tan 8 case, for simplicity we ignore less 
important terms and keep only the leading part contributions given in the following: 


e Charged Higgs (Fig. fi) 


+ + mem tan? B 
CE (mw) = -CE (mw) = = 1 P, (xg) (3.2) 
Ww 


'The most general Hamiltonian in low-energy supersymmetry also contains the operators O}, Q which are flipped chirality 
partners of O;, Q;. However, they give negligible contributions and thus not considered in the final discussion of physical 
quantities 


where zij = m; /mj and the one-loop integral functions P; are given in the appendix. 
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FIG. 2. The Feynman diagrams involving chargino which give dominate contribution to CQ, » 


e Chargino (Fig. B 


X E: n? 
C$, (mw) = B x E (i,k, k’) Da a cos B m 
4,4 —1 k,k'— 
m. 
—óxy T3( ‘i)r4 P3 (ez, iX dar arse ) Ôkk' RS 
1 2,442 
mw 5 — 38in^O0w ; 
me (AES = T3(k, k ) P (vg, ie UE iX; =) 
+k —3 Vin Uia P (2s. myt Tet gt URS y )} ; 
xi 
xt mem V2 mar 
C$, (mw) = 7 n?’ 8 5 D T1 (i, k, k')Uv2 [ooo cem 


28, VA Uis Ps ( 


TE E as 


—Ókk! BASE m 


where a is the mixing angle of neutral components of the two Higgs doublets, r} = cos? omg + sin? om, 
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(3.3) 


(3.4) 


TQ — 


sin 2o(mg? — my?) A = m,/V2mw sing and àp = my/V2mw cos B are the Yukawa couplings of top and bottom 
quarks, respectively. U, V and N are the matrices which diagonalise the chargino and neutralino mass matrices. 
m = mg = mi is defined as the common mass of the first two generation squarks and all sleptons which are assumed 


to be degenerate, Nj = i tan Ow Nji — 
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FIG. 3. The Feynman diagrams involving neutralino which give dominate contribution to Ca, , 


e Neutralino (Fig. B 


'The formula for the contribution from the neutralino are given by 


4 3 
<0 mem, 1 Togo . 
C$, (mw) = 1 Ki Kw tan? B 1 J NITA [s Okk! ms. T4(j, k)r4 Pı (zsoi,. .) 


jj'=1 k,k'=1 


TÓkk; (2r. Tz (55,3, 3) m ro (j^, 3,4) TRP CNN 


TI. o 
Xx; 


+ôkk' = (2r. T5, 9’, 3) — rors (j, 9’, 4)) TB, CN 
xS 4 E I 


mw Fo 7 » 
S M k J)PS CM an) 


el 


(3.5) 


ny Dos 
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e Gluino (Fig. 
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It is noticeable that the supersymmetric contributions to CQ, , have an overall enhancement factor tan? B. Moreover, 


s ; ; dy 16g? 5 : ; ; 
the gluino loop contributions have an additional enhancement factor E Rm So sizable contributions from 


neutral Higgs penguin diagrams are expected for a sufficiently large tan £. 


IV. B RARE DILEPTONIC DECAYS IN THE MSSM 
A. Inclusive decay B — X447 


Neglecting the strange quark mass, and with p standing for the momentum transfer, the effective Hamiltonian B 
leads to the following matrix element for the inclusive b — s€t@~ decay, 


Qem Gr 
2/2n 


+05! sy, = ys )blyyl + Ciosy (1 AN ys )blyu ysl 
*CQ, 8(1 + 75)bE + Co, 8(1 + ys)blys£} . (4.1) 


M= Kaki, f- vp cu " sity yr (1 + y5)b 


The Wilson coefficients can be evaluated from mw down to the lower scale of about mẹ by using the renormalization 
group equation. When evolving down to b quark scale, the operators O; (j = 1— 6), Qs can mix with Oj, (i = 7,9); 
however, they can be included in an “effective” O7,9 because of their same structures contributing to the b > s£*£^ 
matrix element. As for long-distance contribution from the intermediate J/W family, we follow Ref. and include 
the effect in "effective" d f, Expanding C; in powers of as, i.e., C = C? + 2sC1, one obtains the leading order 
effective Wilson coefficients 


8 
€ È 8 — se Qi 
Co^ (mi) = 99/2 (Op (mw) + 32 — 10g (mw) - 0.012Co,(mw) + Dhan, (a2) 


i=l 


E [cd 2 
COT my) = CD (my) + z BCs + Ca +305 + Ce] 


1 1 
7359 (1, s) [4C3 + AC4 + 3C5 + Cs] — 59 (0, s) [C3 4- 3C4] 


m2 30 my, I (V; > H) 
Tg m2) wa" m2, — p? — imy;T v, 
b em V;—w sse Vi i i 

x [3C1 + C2 + 3C3 + C4 + 3€5 + Co], (4.3) 


Cio (ms) = Cio(mw), 
Ca: (mp) = y PS (my). 


Here s = p”/m? is the scaled dilepton invariant mass square, , n = as(mw)/as(mp), vector hi, ai, f; and g; are given 


in [24] and 


Mb 
Cos (mw) = ~~ (Ca, (mw) + Co, (mw), (4.6) 
Ci.g = (—0.4561, 1.0208, —0.0041, —0.0603, 0.0028, 0.0037). (4.7) 


At NLO level, the Wilson coefficient oe f can be found in Zi and 


. 37208 e : 8 : 
C} TF (my) Ni 776/23 (n = Do? Hf (my) apum T my) de arm u gi yc. ff my) 
4 14/23 2/23 1791104 1674721 25/23 | Oef f 
d 421 441 diu sacl , 
pega eee 25 25 "7 |S 
8 
+> [emei (mw) + fi + gin] i (4.8) 
i=1 
e 1 e ; 2 5+4 
Cy fI (gy) = — O) TF (my) | —4Lie(s) — 21n(s) In(1— s) - £r? — : In(1 — s) 
12 3 1-4 2s 
2s(1 + s)(1— 2s) 5+ 9s — 6s? 
—————————— 1 ———————_- | . 4.9 
(1—sPü43s) 9) + 3n — 0425) Ven 
Function g(m2/m, s) in Eq. E3) arises from the one-loop matrix elements of the four-quark operators, and 
4 8 16x 4, 2z 4r V? f In Z(z,y) im. for 4v/y <1 
g(z,y) = -gln Ptr gut Mb darcian a for 4z/y>1 (4.10) 
where 
lega enm 
Z(xz,y) = ————. (4.11) 
1—,/1— B 


parameter & as 2.3 [25] in numerical calculations. 


'To estimate the Yi s pp contribution in the second term in brace of Eq. (1.3). we take the phenomenological 
The formula of invariant dilepton mass distribution has been derived in [13], which is given by 
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Ar 
ds = 768n5 em ) 


Uc KT PaA- s) (1 -— p 1/2 p(s) (4.12) 


with 
ef f 2 2r 2 ef f 2 2r 2 or 
D(s) = 4|C7!? | (1+ 0 2) lcs | (1+ —)(1 + 2s) + [Cio] (1— 8r + 2s+—) 


€ e * 2 3 3 " 
412Re(C8ff pef f*( + =) + 5|Ca, ?(s — 4r) + 5lCo; [^5 + 6Re(CyoCh, )r!? (4.13) 


where r = m2/m. To get rid of large uncertainties due to m? and CKM elements in Eq. (1.13). we normalized the 
decay rate to the semileptonic decay rate 


G2 5 " " 
I(B > Xelv) = TEP Kaf F(T) 


1927 m" (4.14) 


Here f(x) = 1 — 8x? + 8x? — a? — 24rtln v is the phase-space factor and k(x) is a sizable next-to-leading QCD 
correction to the semileptonic decay pe]. 

The angular information and the forward-backward (FB) asymmetry are also sensitive to the details of the new 
physics. Defining the forward-backward asymmetry as 


J d cos 0(d?T /dsd cos 0) — T d cos (dT /dsd cos 0) 


1 2 0 2 E (4.15) 
Jo d cos0(d?T /dsd cos0) + f, dcos 6(d?T /dsd cos 0) 


Arp(s) = 


where @ is the angle between the momentum of B-meson and £* in the center of mass frame of the dilepton, we obtain 


6(1— 4r/s)1/? 


pay i [207 Cj, + C$ Chs + 202 C2, 7? + Off o n2]. (4.16) 
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B. Exclusive decay B, — 44y 


Now let us turn to the rare radiative decay B, — £* (^^y. The exclusive decay can be obtained from the inclusive 
decay b — s£* ^», and further, from b — s£*£^. To achieve this, it is necessary to attach photon to any charged 
internal and external lines in the Feynman diagrams of b — slt. As pointed out in Ref. [3], contributions coming 
from the attachment of photon to any charged internal line are strongly suppressed and we can neglect them safely. 
However, since the mass of (-lepton is not much smaller than that of B,-meson, in B, — £* ^^y decay, the contributions 
of the diagrams with photon radiating from final leptons are comparable with those from initial quarks. When a photon 
is attached to the initial quark lines, the corresponding matrix element for the B — £* (^^ decay can be written as 


3/2 
Gem Gr * * : * * f 
M3 - Vas KiKi { [AE paso e paqo + iB(es (pa) — (“pau ) Cyl 
+ (CE nabo APB + iD(e (pq) m (e*»)au) by yst) , (4.17) 
where 
1 € € 
A=—, lc H Gp?) 2c! Ts?) 
mpg 


D = —-Hh(p). (4.18) 


In obtaining Eq. (4.17) we have used 
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(y|Stop yp (1 - *ys)b| Bs) = m [6noge Podo Ga (p?) i [(e (pq) > (€*p)qu)| F2(p?) } , (4.20) 


and 


(78. + 75)|Bs) = 0. (4.21) 


Eq. ) can be obtained by multiplying p,, in both sides of Eq. (4.20) and using the equations of motion. From Eq. 
(4.21) one can see that the neutral scalars do not contribute to the matrix element Mı. 
When a photon is radiated from the final ¢-leptons, the situation is different. Using the expressions 


Here e,, and qq are the four vector polarization and momentum of photon, respectively; G;, F; are form factors PIBS). 
fal 


(0|50|Bs) = 0, 
(0|so,,, (1 + y5)b| Bs) = 0, 
(Olsyuys|Bs} = ~i fB. Pau (4.22) 


and the conservation of the vector current, one finds that only the operators O4.» and Og give contribution to this 
Bremsstrahlung part. The corresponding matrix is given by p 


3/2 2 
Qem GF ; MB IE Pp, PB ‘| 
Mz = KyKii2 Coe aia e 
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1 l ks = 

BCRA ame + E e+ MEALEN — ty \ (4.23) 
2mem 2p1q — 2paq 2piq 2p2q 


Here Pg., fp, are the momentum and the decay constant of the B, meson, p1, po are momenta of the final ¢-leptons. 
Finally, the total matrix element for the B, — £*£^^ decay is obtained as a sum of the Mı and Mg. After 
summing over the spins of the -leptons and polarization of the photon, we get the square of the matrix element as 
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It is obvious that the quantity |.M|? depends only on the scalar products of the momenta of the external particles. 
In this paper, we follow Ref. and consider the photon in B, — £^£^^y as a hard photon and impose a cut on the 
photon energy Ey Hl which correspond to the radiated photon can be detected in the experiments. This cut requires 
E, > 6 mg,/2 with 6 = 0.02. 

After integrating over the phase space and the lepton energy Ei, we express the decay rate as 


?When photon is soft, both processes of B, > €+@~y and B, — £*£^ should be considered together, and in this case, the 
infrared singular terms in |M2|? can be canceled exactly by the O(@em) virtual correction in B, — £*£^ B. 


3/2 5 2 1—6 
eG cz pb, [hs f "Eae CPC 
= K Kis s E 1— S8)?d8 Dens A +/B §+ 2f 
Van enl ss i fy e Sem [AP + |BI*) + 27) 


+ (ICI? + |D?) (8 — 45)] - 2f. 


1-6 
MCR y ef (1 — 8)?d8Re(A) In 2 


2m 
C 1-6 T 

me DAS] ( (1— 8)d$Re(B) | (1 +4 +58)In 2+ &41— Z 

2Memp $ 
4 2 m? C 1-8 AP — 28 Af. 
-Be org FE pg f dé |(1+8+——)In 2+1- — 
MB, 2mm, Af p E 

C Rd 2 — 10f + 8f? 2(1— 4f)8 Af 
+f3,| pf alasi c ECTS, t ILU 1-2 (4.28) 

4f — Ss — 


where $ = p?/m$, , f = mz /m3., £ = Z(f, 8) takes the form given in Eq. (4.11). 


V. EXPERIMENTAL CONSTRAINTS ON THE RELEVANT SUSY PARAMETERS 


Before scanning the relevant parameter space of MSSM, we assume (i) the masses of the particles are restricted to 
the sub- TeV regime and larger than the lower experimental bounds XE (ii) A-parameters are smaller than 3M6 P9); 
(iii) Except for the third generation squarks, all sfermions are degenerate and have masses of ~ 1TeV; (iv) The GUT 
mass relation Mı ~ M3/2 for gauginoes is used; (v) Only small flavor violation e; < 0.1 is allowed. We also take into 
account the well-known large radiative corrections to neutral Higgs masses [po]. With above assumptions and gluino 
mass fixed as its lower experimental bound 190 GeV P1) , the relevant parameter space is determined by nine input 
parameters: Mo, M;,, M;., p, tan B, At, Ap, Mz and ma. In addition, we consider the following the experimental 
constraints in our scan: 

(1) The recently reported value of muon g — 2 iE] shows a 2.6 standard deviation from its SM prediction. The SUSY 
explanation of this deviation requires (i) u > 0 and (ii) large tan 8. In our calculation we assume 20 < tang < 50 
and at least one of the charginos or neutralinos must be lighter than 500 GeV (ig). 

(2) Non-observation of any supersymmetric signals at CERN ete~ collider LEP-II and the Fermilab Tevatron 
imposes lower bounds as 


My+ > 78.6 GeV mpo > 88.3 GeV, ma > 67.7 GeV, 
mo > 42.0 GeV mz, 2 86.4 GeV, mj, 2 75.0 GeV. 


(5.1) 


(3) The latest measurement of the inclusive branching ratio from CLEO and BELLE B gives world average value 
2.44 x 1074 < Br(B — X,y) < 4.02 x 10^ (95% C. L.), (5.2) 


which is specially useful to constrain extensions of the S T) eun studies used leading-order (LO) SM result to limit 
the MSSM parameter space. Instead of LO calculation [|l , the branching ratio of B(B — X,7) has been estimated 
at next-leading-order (NLO) level in the SM p4 with Ce 22. 596 increase of its central value and theoretical error 
less than 1096. In this paper, we use the SM result computed at NLO level whereas the additional supersymmetric 
contributions at LO level. We do not use the available NLO matching conditions for the supersymmetric particles 
since they are computed under the specific assumptions about the sparticle spectrum, not necessarily satisfied in the 
criteria, and moreover, they are not valid for large values of tan 8 B. Since the Wilson coefficient Cg accounts for 
only 3% of the standard model b — sy amplitude, it is therefore not expected to be significantly more important in 
the MSSM. In this case, Eq. (p.3) implies the constraints on the ratio Rz = Cae. 


0.83 < R; < 1.18, or —1.24< R7 < —0.94. (5.3) 


(4) The stringent bounds on the magnitude of the short distance coefficients come from the Collider Detector at 
Fermilab (CDF) [R1] 
Br(B* > Kt£*€-) < (60,5.2) x 107° (90%C.L.), 
Br(B° — K*9£* £7) < (290, 4.0) x 107° (90%C.L.), 
Br(B, > 4£'£ ) < (5.4,2.0) x 107° (95%C.L.) (5.4) 


for l = e, u. 
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VI. NUMERICAL RESULTS AND DISCUSSIONS 


We perform a complete scan over the nine-dimensional parameter space of the MSSM. For reference, we present 


our SM predictions 


Br(B — X,U £7) = (11.8, 9.87, 2.94) x 107°, 
Br(B, — £* £^) = (2.50, 2.62, 5.49) x 1078 


(6.2) 


for l = e, 4,7. These values are obtained by taking the QCD coupling constant a;(m»,) = 0.218, the masses, decay 
widths and branching ratios of J/W family in Ref. BI, the normalized factor Br(B — Xlv) = 10.2% and the 


form-factors 


1 GeV 3.74 GeV 
G (pP = ————s. QP = 
i) (1— p?/5.62)2' xp) (1— p?/40.5)? 
0.8 GeV 0.68 GeV 
F ee nci DNUS F: 2 = 
1(p ) (1 = p^ /6.52)2 , »(p ) (1 — p?/30)2 


as well as the fixed input parameters listed in Table I in the numerical calculations. 


Table I. The value of the input parameters used in the numerical calculations (mass and decay constant in unit GeV). 


Mt Me Th Me mB, mw 
176 1.4 4.8 1.78 5.26 80.45 
fs, |KoKi| |KoKi/Kol o, T(Bs) sin“ Ow 
0.21 0.045 0.95 137 1.64x107-°**s 0.233 
20 ~ 104. 
Xu (a) 
4^ Tor 
o 7 or 
>s tol 
E M 
o 2 
[aa] 10* — 
f: 8 
a 
m s E 
ra 
j 
10 : 
10 — 
9 E 
8 E 
7 E . . : 
ml à, balai wer] ah al pea 
10° 10 


D 
£1 


(6.3) 


FIG. 4. The parameter space scatter plot of Br(B — X,£* £^) (left) and Br(B, — £*£ y) (right) vs ex for £ = 7 in (a) and 


for £ = u in (b). The solid lines stand for the SM predictions 
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The results shown in Fig. H indicates that in some part of parameter space, the enhancement factors R of the branching 
ratios can be 5 for the inclusive decay B — X,7*7~ and a couple of orders of magnitude over the standard model 
predictions for the exclusive decay B, — r* ^y. This is quite different from the previous studies (1. As illustrative 
examples, the branching ratios for ditau final state as functions of the gluino mass are plotted in Fig. H The effects 
drop with the increase of gluino mass, showing the decoupling property of the MSSM. Table 2 shows the case in which 
both enhancement factors reach their maximum values and the corresponding SUSY parameters. A comparison of 
the enhancement factors with and without the NCL effects is also presented. 


Table II. The maximum enhancement factors and their corresponding SUSY parameters (mass in unit GeV). 


Rmar(B > X,rT) Rmar(Bs > TTY) 
With NCL without NCL] With NCL|without NCL 
m 
tang —40 ma = 453 m; —143 m; = 92 


We have the following comments on the results: 

(1) The constraint on B, — £*£- is of special useful to limit the effects of the scalars and pseudo-scalars. However, 
since no helicity suppression compared with Cio (see Eq. 1.23), the contribution of Cg,, in B —> X,r*r- and 
B, — £* U^» decay is still dominate for large tan £. 

(2) As helicity suppression, the effects of Co, , in B — X,r*r- are very different from those with other dilepton 
final states (see Eq. 1.13). The smaller mass of lepton in the final states, the less effect of gluino and the neutralinos. 
In the inclusive decays with lighter dilepton final states, the dominate contribution comes from the interaction term 


2 Re(C£! j ce f ) with opposite sign to its SM value. For large tan 8 region, relative large m, make the neutral Higgs 
with NCL contribute sizably to the inclusive decay. 


- 
[I 


Branching ratio (10°) 
— 


1 l D 
100 400 00 1000 


Mass of gluino (GeV) 


FIG. 5. The branching ratios vs gluino mass with tan 8 = 40, ma = 271 GeV, m;, = 320 GeV, mz, = 265 GeV. The solid 
and dashed lines correspond to the inclusive and exclusive decays with £ = T, respectively. 


Some numerical examples in the region of R7 > 0 are presented in Figs. B8 Fig. id and Fig. ai show the branching 
ratios of B — X,r^r- and B, — r*7 ^, respectively. For the specified parameter values in the figures, one sees 
that the contributions from gluino and neutralinos are dominant. When the mass of the CP-odd Higgs is less than 
700 GeV or tan B > 40, the branching ratio of B —^ X,r^7- can be enhanced by a factor 2, and for B, — r^r y, by 
one order over the SM results. 
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~ 
© 
o (a) 
w 7L 
~ 
X 6 + 
BPa 
x 

5 L 
= 
tü 

4 BE 

| | 
100 400 700 1000 
m, (GeV) 


FIG. 6. Branching ratio of B — X,r*r- versus ma in (a) and versus tan in (b). 
for the SM prediction, whereas the solid (dashed) one is the MSSM prediction with (without) the contributions of 
The parameter space in (a) is specified as m; = (179,304) GeV,m; = (257,700) GeV, 
Mgt = (94,417) GeV, mgo = (81,109, 200, 417) GeV and tan G = 20; In (b), ma = 899 GeV, the masses of the sparticles are 


the gluino and neutralinos. 


dependent on the value of tan f. 


= 
o 
W 


Br(B c) (10°) 


1 
100 


The differential branchin 
squared are plotted in Fig. 
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m, (GeV) 


Br(B—X,:*c) (10°) 


Br(B , c) (10°) 


- 
Ss 


EN 
e 
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N 


1 
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FIG. 7. The same as Fig. (6, but for B, => T* T». 


dilepton mass squared is shown in Fig. 
predictions, especially in large invariant dilepton mass region. 


ratios of B > X,r^T- 
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and B, — r*T- * versus the scaled invariant dilepton mass 
, while the dependence of the FB asymmetry for B — X,7+7~ on the scaled invariant 
ig. The figures show significant differences between the SM and the MSSM 


The dotted line stands 


dr(B5X;*cyds (10°) 
| PoP rT en 


FIG. 8. The same as Fig. Ba 


We should point out that since some common contributions appearing in both the numerator and the denominator 
cancel out to some extent, the FB asymmetry is a sensitive, relatively model-independent probe of these models. 
We stress that all these distributions would be useful for fitting the future experimental results in the framework of 
the MSSM, especially when some deviations from the SM predictions are discovered in future experiments. Different 
models, such as the MSSM and 2HDM’s, may all predict some enhancements, but they may give different behaviors 


for some distributions. 


0.5 


-0.5 


FIG. 9. The same as Fig. (op 
mass squared s with ma = 305 GeV. 


Although decay modes with the ditau final states are experimentally difficult compared to their di-muonic counter- 
parts, the considered decay modes are more sensitive to new physics. More theoretical studies, such as higher order 
effects which might also lead to important modifications as in the case of B — X, BJ, and experimental efforts for 


the decays considered are valuable. 


dr (B,«*cJ/ds (10°) 


, but for the differential branching ratios of B — X,rT^r and Bs T 
invariant dilepton mass squared with ma = 305 GeV. 


+ 
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, but for the forward-backward asymmetry of B —5 X,r^T versus the scaled invariant dilepton 


T ^y versus the scaled 


VII. CONCLUSION 


In this work we performed a complete one-loop calculation of the inclusive decay B — X,€*@ and exclusive decay 
B, — £* ty in the MSSM. Various experimental constraints on the relevant SUSY parameters, such as B — X.7, 
B —> KOLO, B, —> £*07 and the latest gu — 2 experimental measurement, were considered to constrain the 
parameter space of the MSSM. Our results showed that the contributions from the gluino and neutralino loops, which 
were neglected in previous studies, might be quite important or even dominant in some part of parameter space. 
'These supersymmetric contributions could significantly enhance the branching ratios over the SM predictions. Also, 
with these contributions the distributions of the forward-backward asymmetry of B — X,7*+7~ and some other 
distributions could differ significantly from their SM predictions. Up to now the currently running B-factories such 
as BaBar at ee and BELLE at KEKB B4 have collected about 3.2 x 10° BB pair with a luminosity of 
(3 ~ 4) x 1033 cm7 =t p3lp4. And the BaBar will take 105 BB pur in three years while the BELLE 10’ BB 
pairs each year. Since the branching ratios of the inclusive decay B — X, T+T and exclusive decay B, —^ r* ry can 
be enhanced by the SUSY effects with large tang to reach the level of 1075, these ditau decays might be observable 
at the B factories. 


ACKNOWLEDGMENT 


We would like thank Prof. A. Ali and C.-P. Yuan for useful discussion. This work is supported in part by a grant 
of Chinese Academy of Science for Outstanding Young Scholars. 


APPENDIX 


In the Appendix we present the operator basis O;, Q; and one-loop integral functions for b — s£* £^. 
The operator basis O; is the same as the one used for the b — s£*£^ in the SM, 


0ı = Say" Leg) (Cay Lba), 
O2 = (5a Lea) (Esy" Lbs), 


O35 = Sy" Lew) (= itis) ; 


q 


O4,6 = (5a7" Leg) (= soda) ; 
q 


O- = eae (mR + Ms L)ba Fv, 

Os = T 8a0"" (ma R + msL)tasbaGiy, 

Os = Goa Lbs) Paul), 

O10 = T (Say Lba) (C1758), es 


where the chiral structure is specified by the projectors L, R = (1 F y5)/2, while a and 8 are color indices. F „y and 
Giiv denote the QED and QCD field strength tensors, respectively. ta are the color triplet generators, g and gs stand 
for the electromagnetic and strong coupling constants. 

Operators Q; come from exchanging the neutral Higgs bosons in MSSM and are defined by 


Qı = <5 (Ga Rba) (20), 


g? 
Q = Toza Ge Rho )(£s£), 


g? 
Q34 = cs = 5 (SaRba) (Est (R, L) 2] 
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2 
O56 — TE (Sq Rbg) (= LL , 


2 
g v - 
O78 = 1642 Sad” Rba S sg(ft, 2 , 


Qo 110 = T 8,0"" Rbg (= L1 . (7.2) 


q 


'The one-loop integral functions which appear within the MSSM matching conditions are given by 
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P,(z) = — Ins, 
Pio) = cm [zm me = ny) 

Ps(z, y) TEES ET r- my], 
Py(x,y,z) = erences (ao y) + (z« 2] (7.3) 
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